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We explore the superhorizon generation of large /nl of the local form in two field inflation. 
We calculate the two- and three-point observables in a general class of potentials which allow for 
an analytic treatment using the 5N formalism. Motivated by the conservation of the curvature 
perturbation outside the horizon in the adiabatic mode and also by the observed adiabaticity of 
the power spectrum, we follow the evolution of /nl 3 ' until it is driven into the adibatic solution 
by passing through a phase of effectively single field inflation. We find that although large /{j ™ 1 
may be generated during inflation, such non-gaussianities are transitory and will be exponentially 
damped as the cosmological fluctuations approach adiabaticity. 



I. INTRODUCTION 

Amongst the striking successes of cosmic inflation, the 
most dramatic is the natural prediction of the cosmolog- 
ical fluctuations observed in the cosmic microwave back- 
ground and in large scale structure ^E]. Many models 
of inflation predict a nearly scale- invariant, nearly gaus- 
sian power spectrum which accords with current obser- 
vation, though measurements of the scalar spectral index 
and bounds on the relative gravitational wave amplitude 
have begun to rule out some models of inflation [7 . While 
ongoing and future satellite experiments will greatly im- 
prove the precision of measurements of the power spec- 
trum [8l [9] (and thus help to further constrain theories of 
the early universe), there is a limit to the insight gained 
from the power spectrum alone [T01 E] • 

To go further, we can consider the non-gaussian parts 
of the spectrum, which in principle contain a much richer 
set of information, and thus may provide a key obser- 
vational tool to distinguish between detailed models of 
inflation. In particular, there has been a great deal 
of progress in measuring the bispectrum (3-point func- 
tion) from the CMB [7] and from large scale structure 
|12j . Complete characterization of the 3-point function 
is a difficult task, but a useful basis can be given in 
terms of different shapes of momentum-space triangles. 
In particular, one considers the following forms (we de- 
fine ki > k-2 > fca in the following): local/squeezed 
(fci ~ &2 ^ kz), equilateral (k± ~ &2 ~ fe) and orthogo- 
nal (constructed to be nearly orthogonal to the preceding 
two forms) - see [T3] for a review of the theoretical and 
observational motivation for this basis. Of course, these 
are not the only forms one can write down, and other 
shapes may probe other aspects of inflation. 

The current observational bounds are — 5 < /^L al < 59 
(WMAP7+SDSS), -214 < /°f 11 < 266 (WMAP7) and 
-410 < /N L thog < 6 (WMAP7), where / NL character- 



izes the size of the bispectrum compared with the power 



spectrum (see ( 16 ) for the precise definition) 
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The local shape is especially important because it is 
predicted to be small in all models of single field infla- 
tion [T4Ttl"?] , so a convincing detection of significant local- 
shape bispectrum would therefore rule out a wide class 
of models including the simplest single field inflationary 
scenario. This result does not, however, hold for multi- 
field inflation; in particular there are a number of models 
that use the continued evolution of the curvature pertur- 
bation outside the horizon in order to generate a large 
/^°£ al ; for example see [HH32]- 

This superhorizon evolution of the curvature pertur- 
bation is characteristic of the presence of non-adiabatic 
fluctuations in the cosmological fluid, which is itself one 
of the most important differences between single field in- 
flation and models with multiple dynamical fields. Since, 
regardless of the contents of the universe, there always ex- 
ist two purely adiabatic modes in which the superhorizon 
curvature perturbation is conserved [331 EU j the (only) 
two scalar fluctuation modes in single field inflation are 
guaranteed to be adiabatic and remain so for the subse- 
quent evolution of the universe |35j . This means that for 
a given model of single field inflation, we need only cal- 
culate the evolution of the perturbations until the time 
of horizon exit in order to make predictions about the 
observational effects of cosmological fluctuations. Mod- 
els with multiple dynamical fields, however, admit non- 
adiabatic solutions in which the curvature perturbation 
can evolve outside the horizon. As a result, we must cal- 
culate the evolution of the curvature perturbation until 
it becomes conserved, or until it is observed. If the cur- 
vature perturbation does not become conserved, models 
with multiple dynamical fields cannot make sharp pre- 
dictions about cosmological fluctuations unless we have 
a complete understanding of the entire history of the uni- 
verse stretching back to the time of inflation. 

If there are non-adiabatic fluctuations present in the 
cosmological fluid during the radiation-dominated era 
this would leave observable effects on the cosmic mi- 
crowave background. Specifically, these fluctuations 
carry a different phase than adiabatic fluctuations, and 
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have peaks at different positions in the angular power 
spectrum [36 39J . Measurements of the power spectrum 
can place upper limits on the magnitude of the non- 
adiabatic contribution to the spectrum of fluctuations 
in the early universe, and current observations have not 
found any evidence for such a contribution [7]. 

The non-observation of a non-adiabatic (or isocurva- 
ture) contribution to the spectrum is not by itself enough 
to rule out multifield inflation, since there are at least 
two ways in which the perturbations produced by multi- 
field inflation can become adiabatic. The first is to pass 
through a phase of effectively single field inflation, by 
which we mean that all but one of the eigenvalues of the 
mass matrix for the inflaton fields are large and posi- 
tive. Secondly, the adiabatic solution is attractive if the 
universe passes through a sufficiently long period of local 
thermal equilibrium with no non-zero conserved quantum 
numbers [301 I3T] . 

Now, as observed above, the principal suggestions for 
having a large /nl 3,1 generated in a multifield model use 
the evolution of the curvature perturbation outside the 
horizon in the presence of non-adiabatic fluctuations to 
boost a initially small bispectrum. In this way, super- 
horizon evolution of the curvature perturbation allows 
for the generation of non-gaussianity, even if fluctuations 
were nearly gaussian at horizon exit. 

While the above scenario is apparently promising with 
regards to the generation of large f^ 31 , ^ * s a l so m " 
complete. If there remain non-adiabatic fluctuations in 
the cosmological fluid, the curvature perturbation and 
its correlations (including the bispectrum) continue to 
evolve. Thus any calculated non-gaussianity in a partic- 
ular multifield model is completely unpredictive unless it 
can be carried through to an epoch in which the pertur- 
bations either become purely adiabatic (i.e. a phase of 
effectively single field inflation or a period of local ther- 
mal equilibrium) or are observed. 

Our goal in this paper, then, is to study the contin- 
ued evolution of f^£ al until the curvature perturbations 
becomes constant and thus address whether there exist 
models of multiple field inflation which predict observ- 
ably large local non-gaussianity and a purely adiabatic 
power spectrum. We insist on the latter not only to be 
consistent with observational bounds, but also to be cer- 
tain that any alleged non-gaussianity would survive to 
the present epoch. 

To this end we use the SN formalism to obtain ana- 
lytical formulas for two- and three-point statistics in a 
wide class of potentials for two-field inflation. We then 
require that the fluctuations become purely adiabatic be- 
fore inflation ends by passing through a phase of effec- 
tively single-field inflation. For the potentials we study, 
the suppression of non-adiabatic fluctuations by a phase 
of single field inflation also dynamically suppresses the 
local bispectrum. We therefore conclude that for this set 
of models, regardless of initial conditions, a purely adi- 
abatic power spectrum and large local non-gaussianity 
cannot be simultaneously produced from inflation alone. 



Our conclusions about /^l 3,1 seem to be in agreement 
with numerical studies of particular models of multifield 
inflation; see for example [55J 03] . 

Of course, the space of models that could produce non- 
gaussianities and an adiabatic power spectrum is consid- 
erably larger than the set we examine. In particular there 
are models such as the curvaton scenario and modulated 
reheating, where the non-adiabatic fluctuations persist 
through the end of inflation and are damped away during 
a period of local thermal equilibrium. It is reasonable to 
suggest that in these models, or in a different class of po- 
tentials, non-gaussianities generated from superhorizon 
evolution of the curvature perturbation may be observ- 
able. We will discuss these possibilities below and will be 
exploring them in more detail in subsequent work. How- 
ever, we believe that without understanding the complete 
evolution of perturbations until they become adiabatic, 
it is somewhat premature to claim observable as a 

prediction of multifield inflation. 

The structure of the paper is as follows. In section [n] 
we review the SN formalism, find the form of the most 
general two dimensional potential to which such formal- 
ism can be readily applied and derive expressions for the 
statistical properties of the curvature fluctuations for two 
particular forms of our potential (the painful details of 
the general case are relegated to appendix |B| . In section 
III for the same two potential forms (again the general 
case can be found in appendix [b]), we show that a phase 



of single field inflation drives f^f^ 1 to be slow roll sup- 
pressed. In section |IV| we discuss how our results may 
be circumvented. And, finally, we summarize our conclu- 
sions in section W\ 



II. THE MODEL 



A. SN and the Spectrum 



Consider an action of the form: 



S 



d 4 x 



\ m l R 



(i) 

We concern ourselves with two field models with canoni- 
cal kinetic terms, though, as we note below, some exten- 
sions to more fields will not substantially change our cen- 
tral conclusions. In this case G a b = S a t, with a,b = 1,2, 
and for clarity we use <p and \ for the two scalar fields. 
The equations of motion are then given by: 



H 2 



= + 3Hcj> - 
= X + 3Hx 

1 n -4 

2 ' 



Zrrvi 



-d x W 

+ \x 2 + w {<t>,x) 



H 



2m 2 p 



r + r 



(2) 
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Assuming slow roll (\H\ <C H 2 and \(f> a \ <C H\<f> a \) these 
become: 



3if^ 
H 2 



-d x w 



3m 2 



W . 



(3) 



Models of this type and their associated non- 
gaussianities have been well studied in the literature, for 
specific classes of potentials - for example sum-separable 
potentials (e.g. 02), product-separable potentials (e.g. 
gH HI]) and many others (e.g. [22H29]). Such stud- 
ies mostly use the SN formalism [4"6K4T)] to calculate the 
evolution of the curvature perturbation, £. Accordingly, 
we begin by reviewing this formalism - our treatment 
follows closely that of Vernizzi and Wands [42] . 

C is defined as the curvature perturbation on uniform 
density hypersurfaces [50] and is of particular utility, 
since it is conserved for adiabatic modes (i.e. those from 
single field inflation) whose physical size is greater than 
that of the horizon. In terms of the perturbed metric 
(tensor and vector perturbations are suppressed): 



ds 2 = -(1 + 2A)dt 2 + 2B ti dtdx 
+ a{t) 2 [(I - 2iP) 5 t] 
C is given up to second order by: 



2E. ij ] dx l dx 3 



(4) 



C = -ip-ip 



H r 1 ; r H r r" 1 

— op+-tpdp+—dpdp + 



H 

2p [j 



5 P \ 
(5) 

where p is the energy density and 8p is its perturbation. 

The heart of the SN formalism is the observation that, 
at large scales, £ (at some comoving time t c ) is given 
by the perturbation to the number of e-foldings from an 
initially flat hypersurface (at t = t*) to a comoving one 
(at t = t c ): 

((t c ,x)~6N(t c ,U,x)=Af(t c ,t*,x)-N(t c ,t*) . (6) 

N is the unperturbed number of e-foldings, given by in- 
tegrating H from t* to t c : 



N 



Hdt. 



(7) 



If we take t* as the time a particular mode exits the 
horizon (i.e. when k = aH), the number of e-foldings 
can be viewed as a function of the field configuration 
on the hypersurface defined by horizon exit, (f) 1 (t* , x) , 
and t c . The perturbation in TV can then be expressed 
(up to some order, in our case second) in terms of the 
fluctuations of the scalar fields at horizon exit: 



5N 



E 

i 



N,iS<t>i + Nut 



(8) 



i j 



The derivatives here are with respect to the fields at t = 
t* (N,i = §^r). In reality, SN should also depend on the 
values of the scalar field velocities. However, so long as 
the slow roll approximation holds at horizon exit, only 
the field values determine the subsequent dynamics. 



Statistics 

The utility of the SN formalism is that it provides a 
readily tractable way to calculate cosmological observ- 
ables; in particular, two and three-point statistics. We 
summarize the relevant results below, again we closely 
follow the treatment of Vernizzi and Wands [12] (see also 

The power spectrum, V( is defined by: 

2tt 2 

(c fel a 2 ) = ^^(fcx + fe 2 )^p c (fc 1 ) . (9) 

In a similar fashion the two-point correlation function for 
the scalar field fluctuations at horizon exit is given by: 



2ir 2 



p*{k) = 



Ml 

4^2 



Then, from (g, @ and ([To]), we have: 



(10) 



(11) 



The spectral index is given by (the approximate equality 
denotes lowest order in slow roll): 



l _ dlnP c _ dln"P c _ 1 



dln7> 



dlnfc 



AN 



H dt 



Using the above definition and ( 11 ) we have 



»< 



1 = -2e 



2 EuhNuNj 



(12) 



(13) 



e = —H/H 2 is the usual slow roll parameter. Using the 
slow roll equations of motion, this can be written as [52j : 



-2e 



2m 2 p J2 I jWjjN J N.j 



WJ2 k N 



K ly2 K 



(14) 

Three point statistics can be obtained in much the 
same fashion. The curvature bispectrum, B^, is defined 
through: 

(CfeiCfe.Cfea) = (2tt) 3 S^ (j2 k ^J B ( (kuk 2 ,k 3 ). (15) 

The bispectrum can be used to define the non-linearity 
parameter /nl (defined in 14 ) which is the quantity 
most often referenced in observational constraints: 



5 Eikfin^ 



(16) 
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From ^ and @: 

(CfexCfesCfes) = E N ,i N ,jNk WkMkJ^t) 



UK 



IJKL 



perms. (17) 



The star denotes a convolution and higher order terms 
have been neglected [5TJ [53] . 

The three-point field correlation for the fields is given 
by |5U: 



\ z / j. xi i perms 



(18) 



The sum is over simultaneous rearrangements of /, J and 
K and the momenta k±, k 2 and /C3 and M. is given by: 



M(k 1 ,k 2 ,k 3 ) 



h 2 k 2 1 ^2 1,2 



(19) 



Carryin g ou r the sum over permutations gives the first 
term in (jlTj) [22] as (using J2i ^.i4>i = H): 



Y,N,iNjN k Y <i>i5jKM(k u k 2 ,k 3 ) = 

Hj2 N }J r (ki,k 2 ,k 3 ), (20) 



UK 



perms 



The function / is given by: 

f(k 1 ,k 2 ,k 3 ) = -1 



T 



(24) 



and is a shape-dependent (in momentum space) function, 
with values between and 5/6 [H]. Then, defining the 
scalar-tensor ratio r = %V*jn&P^ we can rewrite the 



first term of ( 23 ) as 



6 



(25) 



which constrained by observations to be small since r <C 1 



The second term in ( 23 ) , /ijr , is momentum 



independent and local in real space (so it contributes 
to /nl a ' ) 1 an d can (in principle at least) be larger than 
unity: 

(26) 



6 , (4) = EijNjNjNjj 

5 /NL- {EkN 2 k)2 



The current bounds are -5 < f^l* 1 < 59 
(WMAP7+SDSS), and these should tighten considerably 
with the current generation of CMB experiments, with 
the Planck satellite expected to give A/^f; al ~ 5 [54] . 
As noted above /^L al * s slow r °U suppressed in all mod- 
els of single field inflation and thus could be (in princi- 
ple at least) a sharp probe of the number of dynamical 
scalar degrees of freedom during inflation. However, as 
we shall demonstrate below, generating and preserving a 
large f^ 31 through to adiabaticity is far from easy. 



where: 

F(k 1 ,k 2 ,h) = Y M{k u k 2l k 3 ) 



perms 



= -2 



1 4 £i> ? fc?fcj? -iX>«- 



ki 



(21) 



For the second term in (17), we follow [32] [FI] and 
neglect the connected part of the four-point function and 
use Wick's theorem to obtain (after a few manipulations): 



B(ki,k 2 ,k 3 ) 



4tt 4 P 2 



-1 



C U t kf y4m 2 J2 K N 2 K J2 i kf 
, ErjNiNjNu 



N 2 



From the above and (16), /nl is: 

EijNjNjNjj 



5 /nl 



2m2-p c 



(1 + /) + 



(22) 



(23) 



B. The Potential 

In order for the above expressions to be useful we need 
to take derivatives of N with respect to the initial values 
of the fields. In multifield inflation this is not straight- 
forward, for, unlike in single field inflation, there are an 
infinite number of possible inflationary trajectories. As 
a result, changes in the initial conditions affect not only 
how far along an inflationary trajectory the field has trav- 
eled at some later time, but also the particular trajectory 
the field is actually on. This easiest way to deal with this 
difficulty is to find a constant of the motion that allows 
us to relate the initial and final field values to one another 
- in much the same way as conservation of energy and 
momentum in particle mechanics allow us to connect ini- 
tial and final states without having to know the detailed 
form of a particular solution. 

To construct such a constant we define a vector A (in 
field space) and consider the quantity: 



A- 



(27) 



The integral is defined along a path 7 in field space. In 
order for C 7 to be a good constant of motion it needs to 
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be independent of the particular path 7 (and of course 
constant along trajectories); figure [l] illustrates the con- 
struction of C . In this case we have: 



A-> 



The derivative of C is then given by: 
AC 



At 



(28) 



(29) 



Setting this equal to and using the slow roll equations 
of motion ^ gives: 

(30) 



and thus we can write: 



-d x W x W. 



a.* = 



g{4>,x) 
d x w 



(31) 



For such a \& to exist, we must have (subscripts denote 
derivatives) : 



9x _ ff^W 
Ws Wl 



94> , aWfa 



w 2 



(32) 



In principle any W satisfying the above equation for an 
arbitrary g would allow us to define a suitable C . How- 
ever, if we wish derivatives of C to only depend on the 
choice of trajectory and not the fiducial point at which 
our integral is anchored {((/>o,Xo) in figure [T]) then 
needs to be only a function of <j> and likewise d x ^ needs 
to be only a function of x so that the dependence on end- 
points in C is separable: C — Ci(<p(t),x{t)) + C2(<f>o,Xo)- 
This is turn means both the left and right-hand sides of 
(32 1 are identically zero. 



We can solve the above equations with W(4>, x) = 
F{U{4>) + V(x)) and g = F' . We stress that this is not 
necessarily the only potential for which a suitable con- 
stant of motion can be found, but it is the most general 
one we have been able to construct. With W of this form, 
C is given by: 



C 



1 



1 



V'(x') 



d X ' • (33) 



We note at this juncture that this class of potentials 
has also been recently analyzed in [SS] by Wang (who also 
considers a possible non-canonical kinetic term for the x 
field). Our analysis differs somewhat in the particular 
approach to the calculation. 

With the above form of W the slow roll equations of 
motions can be written as: 



3H<f> ~ F'U' 
3H X ^ F'V' 

H 2 1 



3?7Jp 



(34) 




FIG. 1. The dashed lines show possible classical paths of 
the inflaton (N increases along the path). The quantity C, 
defined as an intergral along the solid lines emanating from 
the point (<j>o, Xo), depends only on the endpoints of the curve 
and does not change as one moves along a trajectory. 



Here, as below, a prime denotes the derivative with 
respect the argument of the function. The slow 
roll parameters (q = {m 2 /2){dW/d(j) i ) 2 and ry.y = 
m 2 (d 2 W/ d(j>i4>j)) are given by (we drop the repeated in- 
dex on the diagonal 77 for conciseness): 



F'U' 



n 



4>x 



F 

F'V 
F 

H 2 

F"U' 2 + F'U' 
F 

F"V' 2 + F'V 
F 

F"U'V'\ „FF" 



F J F' 2 
As noted above, N is given by the integral: 



N 



Hdt. 



(35) 



(36) 



In order to be able to take derivatives of N with respect 
to 4>* and x* (necessary to use the SN formalism) we use 
the slow roll equations of motion to rewrite the above 
expression for N as: 



N 



1 



2m 2 P J* W<t,(4>,x(<l>))' 



2m 2 



wjxAix)) 



A X - (37) 



We've explicitly shown the dependence of each of the 
fields on the other to emphasize that the integrals can 
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only be performed for a subset of all possible trajectories 
- those for which there is a one-to-one mapping from <f> 
to x- Even with this restriction, the calculation of the 
derivatives of N is somewhat involved; for completeness 
we include the gory details in appendix [B] but limit our- 
selves to a simpler (more restricted) class of models in 
the main body of this work. 

The aforementioned simplification comes from further 
constraining the form of the potential, so that N sepa- 
rates into distinct integrals over <f> and x- To do this we 
rewrite (37) as: 



N 



1 



1 F(U(<j>) + V( X )) 



U'(4>) F' (U(4>) + V(x)) r ' 
Then, if we require that F is of a form such that: 
F(U(tf>) + V(x)) 



F>{U{<t>) + V{ X )) 
a little manipulation gives: 

i r i 



u(4>) + v( x ) . 



(38) 



(39) 



N = — 



U'(</>) 



u(4>)&4>- 



v(x) 



V( X )d X 
(40) 

One can show that the most general solution for F of 
this form is given by |56j : 



F(U^) + V(x)) 
F> (U{4>) + V(x)) 



a{U{4>) + V{x)) + W, (41) 



with a, f3 constants. With F constrained thusly, the most 
general form of W is given by: 



W = W a [a{U{4>) + V{x))+2p\ 



l/a 



(42) 



With suitable redefinitions of parameters and potentials 
this can be rewritten as (a ^ 0): 

w = M<p) + v{ x )V . 



So that N is then given by: 

1 r c u(4>) JJL 



N 



7mj 



W (0) 



c v(x) 

V'{x) 



(43) 



d X . (44) 



With the potential in this form our approach is essen- 
tially the same as the one taken by Wang in [55] (see also 
[57]). The case with a = has to be treated separately 
and the potential takes the exponential form: 



W = ^ Exp [U(<j>) + V( X )] 



(45) 



with N given by (the factor of 1/2 comes from the nor- 
malization of the exponent): 



N 



1 



2to^ 



1 



1 



u> WO 

In fact one could just write this as: 



V (x) 



N 



1 



1 



U'((/>) 



d x . (46) 



(47) 



with a similar expression for x- We keep the symmetric 
form for ready comparison with the homogeneous case. 



The End of Inflation 

In addition to the restrictions on the potential (and 
possibly the trajectory) we also want the inflationary pe- 
riod to end with a phase of single field inflation in order 
to damp out the isocurvature perturbations and leave an 
adiabatic spectrum of fluctuations. During such a phase 
the isocurvature perturbations must become heavy while 
the adiabatic ones stay light. The relevant parameters 
are the effective masses of the adiabatic and isocurvature 
perturbations, given respectively by rf a and rf s (where 
a and s are the directions parallel and perpendicular to 
the infiaton motion) [55] : 



e x 7 f _ 2\/lf&Tf x + e^rj x 



(48) 



When studying the fate of /nl in scction [ni] wc will study 
what happens when rj ss becomes large, while rf a stays 
small. 



The Infiaton Trajectory and Violating Slow Roll 

Before moving on to our results, let us, for complete- 
ness, briefly consider what happens if trajectories devi- 
ate from slow roll. To derive the above the results (and 
those that follow below) we have worked in the slow roll 
regime. However, even if we insist on slow roll at horizon 
exit and as we approach the end of inflation, it is easy 
to imagine that slow roll may be violated briefly during 
either a sharp turn in field space (turning being both a 
hallmark of multifield infiaton [35] and necessary for su- 
perhorizon generation of non-gaussianity 19J) or during 
the approach to adiabaticity. To this end, let us consider 
what sort of violations of slow roll may be still tolerated 
in our framework. For a different approach which utilizes 
the Hamilton-Jacobi formalism to study the evolution of 
the Hubble parameter without relying on slow roll, see 
[231 EH [59] . 

To begin with let us consider TV. We stipulate that slow 
roll must be imposed at horizon crossing (t — t*) and 
at comoving surface on which we calculate our spectra 
(t = t c ). Then we consider how the derivatives of TV 
change if we violate slow roll in some region t\ < t < <2, 
where t* < t\ < t2 < t c . Then, writing 



N 



Hdt 



Hdt 



Hdt. 



(49) 



one can show, with a little work and using the fact that 
slow roll conditions apply at t\ and t-i (but not between 
them) : 



dN _ dN 
<9(f>* d(f>* 



SB 



d ( H 



(50) 
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Here, SR denotes the slow roll result (given in appendix 
[A] ) . There is a similar result for ON/ d\* ■ 

For C, we need to check under what circumstances it 
remains a good constant of motion outside the slow roll 
regime. To this end, consider the time derivative of C 
without imposing slow roll. From ( 33 1 : 



dC 



p mu> 

,2 



l Py, 



,F'V' + x 



p 3HV 

,2 X 



m„ m„ 

P 3HU' P 3HV 



(51) 



Thus for our framework to be valid we must satisfy the 
following conditions: 



zero if the friction term in ^ is sub-dominant so that 

4>/U' = x/V = F' . 

To summarize, then, violations of slow roll during the 
inflaton motion are consistent with our SN analysis so 
long as either the violation is sufficiently short (in e- 
folding time) or if H/<j) is only weakly-dependent on the 
initial field value (attractive trajectory) and friction is 
subdominant during any non-slow roll regime. While it 
is unclear how general these scenarios are, we close this 
section by noting that there is some numerical evidence 
to lead us to believe that we can go beyond slow roll. For 
example, Vernizzi and Wands compare analytic and nu- 
merical results for doubly quadratic inflation in [42] (see 
also [43]), for a variety of different initial conditions, and 
they find close agreement. 
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« 1, 



dx< 1, 



mp mv^ x 



(52a) 
(52b) 
(52c) 



The question we should ask, then, is: To what degree 
can we soften the slow roll criteria (\H\ <C H 2 and \(f> a \ <C 
H\(p a \) and still satisfy the criteria of (52 )? If we rewrite 
(52a I and (52b) so the integral is over the number of 



e-foldings, we have: 



C. Results 

We will present our results for both homogeneous and 
exponential potentials here; the involved details of the 
calculation can be found in appendix |Aj but the basic 
procedure is to take derivatives of N and use the deriva- 
tives of C to relate variations of 0* to variations in <p c 
(and likewise for x)- As mentioned above have also stud- 
ied the general potential W(cj>,x) = F(U((j>) + V(x)), the 
details can be found in appendix |B| 



Homogeneous Potential: W(<f>,x) = Wi't 1 ) + ^(x)]" 





d H~ 







diV<C 1 



(53) 



with a similar equation for x- From this we can immedi- 
ately see that any sufficiently short violations of slow roll 
will, as one would expect, leave the equations for deriva- 
tives of ./V unchanged. This is, of course, unsurprising. 
However, the inequality can also be satisfied if the quan- 
tity H/4> has only a weak dependance on compared 
with its magnitude. This latter requirement can be met 
if the trajectory before entering the region of slow roll 
violations is sufficiently attractive, so that <f> and H are 
nearly independent of the initial field value. Unfortu- 
nately, unlike single field inflation, attractor behavior is 
not universal in two field inflation, so we can make no 
precise statements about our model that are potential 
and trajectory inde pend ent. 

With regards to (52c), the most straightforward way 
for violations of slow roll to leave our results unchanged 
is, again, for said violation to be suitable short in dura- 



tion. Alternatively the right hand side of (52c) will be 



Results here are expressed in terms of derivatives of 
the potential (as before a prime denotes a derivative with 
respect to the argument of the function - <fi for U, x f° r 
V and U + V for F) and the slow roll parameters defined 
in (35). First, the mass parameters for the adiabatic and 



isocurvature fluctuations are given by: 



+ 4 (7^1) e e x +£ x„X 



e x v 4> _ 4tLzll £ 4> e x + e <t> v x 



(54) 



For the other quantities of interest, we first define: 
1 



Xh 



Vh 



u m + v, 
1 



V* 



Vc4 ~ U c e} 



Vc4 ~ U c e* 



Then the observables and — 1 are given by: 



(55) 



247r 2 r 



p 



4 

-2e* - - 
7 



1 - 



2et 



- 7 + 1 )X h 



Vh 



1 - 



7 + 1 )Uh 



4 i3 



(56) 



and /nl b y : 



6 (4) _ 2 
5 /nl ~ 7 



7 + 1 W 



if "7 + 1)1* 



(J/.+V.) 2 \ et 4 J *c \ fc 1 

(f + §) 2 



(57) 



Exponential Potential: W(<j>,x) = W o fep[(7(0) + V(x)] 



The slow roll parameters ( 35 ) take a particularly sim 



pie form for an exponential potential of the form of (45 ), 



and readily lead to the following expressions for the mass 
parameters ?7 SS and if": 



n 



(58) 



Next, in a similar fashion to the previous section, we 
define: 



x P — 1 



y e = l - 



(59) 



Then the observables 7\ and tiq — 1 are given by: 



V c = 



96Tt 2 m 



7l£ — 1 



and g/$ by: 




1 e§ Vc \et 4) X e\ \4) 2 J »* \ (e*) 2 ) 



(61) 



III. DAMPING AWAY ISOCURVATURE AND 
THE FATE OF / N l 

Homogeneous Potential: W((j>,x) = [^(<W + ^(x)] 7 

We are now in a position to assess the magnitude of 
fl^ , beginning with the homogeneous potential. Follow- 
ing the arguments of |42j . we see that Xh + Vh = 1 and 



m&x(xh,yh) < 1, where and are defined in (55), so 
the denominator of (57) is of order e~ 2 , where e refers to 
a generic first-order slow roll para met er. Both terms in 



the numerator in the first line of (57) are of order e 



so the first term in parentheses is of order The pa- 



rameter 7 which appears in the denominator of (57) is 



in principle unconstrained, so it might be expected that 
choosing 7 to be smaller than could produce large non- 
gaussianity. However, a factor 7" 1 also appears in (56) 
for the scalar spectral index — 1. Observation of a 
nearly scale-invariant power spect rum [7] thus requires 
that |7 -1 | < 0(1)- Returning to (57), we thus see that 



the second line gives the only contribution to /i7 which 
is not automatically slow roll suppressed: 



6 (4) 



o) 2 ( Xh 

'.) 2 \4 



(*l 



4> e x 



0(e,)+0{l) x 



We are interested in the value of // 4 ^ 



(62) 



NL only after it be- 
while there may large non- 

(4) 



comes a conserved quantity 

gaussianities at any given time during inflation, is a 
dynamic quantity and we must follow its evolution until 
it becomes constant and we can readily translate the pri- 
mordial value into observables today. Alternatively, one 
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could also follow the evolution up to the time of obser- 
vation itself, but that would require detailed knowledge 
of the cosmological history from the present all the way 
back to the inflationary epoch. Given our ignorance of 
the details of much of the early universe, we focus on the 
former approach. Accordingly we wish to study / N l as it 
is forced into a constant value and, since correlations of 
the curvature perturbation will evolve outside the hori- 
zon as long as non-adiabatic fluctuations are present, this 



means we need to evaluate ( 62 ) after any non-adiabaticity 



has been damped away and only the adiabatic mode is 
left. 

As discussed above, we are going to damp away non- 
adiabatic fluctuations by passing through a phase of ef- 
fectively single field inflation, which requires that r] ss be- 
comes large during a finite period before the end of infla- 
tion. Now, in the slow roll approximation the equations 
of motion for the adiabatic and isocurvature fluctuations 
decouple [55], and we have for the latter: 



W 

Ss + 3HSs + -^-f] ss 5s = 0. 
rap 



(63) 



The solution of this equation in the slow roll approxima- 
tion is f58l f60l 



5s oc a(i)- 3/2 



k 
off 



where v is given by 



So for T) ss > I, we find 



9 Wr] ss 
4 m 2 p H 2 ' 



\Ss\ oc a(t)- z ' 2 



(64) 



(65) 



(66) 



and the isocurvature fluctuations are rapidly damped 
away. 

We will now examine the conditions for large i] ss (de- 
fined in (48)). Recall that (35) gives 



J>x = 2 



(7-1) 



(67) 



and I7I must be 0(1) or larger to guarantee scale invari- 
ance. Thus ?]^ x cannot give the dominant contribution 
to ?7 SS when 



n 



e x v 4> _ 4(1=1) e <t> e x + e <P v x 



(68) 



is large. As a result, we must have either ff > \ or 
r] x > j^. If both rf' and r\ x are large, then rf" 7 will also 
be large, and inflation will quickly end before the non- 
adiabatic modes have been damped away, so we will not 
be interested in this case. 



Because rf s appears in the second line of (57), it might 



large non-gaussianity is produced in the process of damp- 
ing away non-adiabatic fluctuations. However, the coef- 
ficient of 77 s3 , specifically ^-§-, is a dynamical quantity 
whose time-dependence is affected by the behavior of rf s . 
To find this time-dependence, we first we note that: 



H 
c'x 
H 



— = -2e V - 2V^x^ x + 4(e^) 2 + 4e^e* (69) 
-2e x r) x - 2^f^rf x + 4(e x ) 2 + 4eV. (70) 



Now, if we have 77* > 1, we can neglect all but the first 
term on the right hand side of (69), and so we find that 



e (t) oc Exp 



-2 / H^dt 



(71) 



Similar remarks apply to e x , in the case that r\ x > 1. 

As a result, we find that after rf s becomes large, the 
dominant contribution to decays exponentially with 
time, with a decay constant greater than that of the 
isocurvature fluctuations. Specifically, from (62) and (66) 
we have for rf s > 1: 



Exp 



Hdt 



2 / C„Hr) ss dt 



(72) 



be expected that for sufficiently large r; ss , observably 



C n is a number which is always greater than 1 whose 
value depends on the particular direction of the effective 
single field during this phase. We thus conclude that 
/nl wu ^ a l wavs be slow roll suppressed upon entering 
the purely adiabatic solution after a phase of effectively 
single field inflation. 

There are at least two cases in which the terms labeled 
C(e*) may be significant, even though those terms are 
proportional to slow roll parameters. In the first case, 
these terms are initially small, but may become large for 
finely tuned trajectories in field space |61j . In this case, 
these terms will always be small after passing through a 
phase of effectively single field inflation for reasons sim- 
ilar to those discussed above. In the second case, the 
individual terms contained in those labeled C(e*) are sig- 
nificant, but there is initially a cancellation so that 
is small at horizon exit [55] (the model studied in [55] has 
a large number of fields, but a similar mechanism should 
work with two fields and an appropriately chosen poten- 
tial). After passing through a phase of effectively single 
field inflation, various terms are damped, and there is no 
longer a cancellation near the end of inflation. This has 
the possibility for giving a significant (though not para- 
metrically large) /nl and purely adiabatic fluctuations 
at the end of inflation. 
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Exponential Potential: W{<j>, x) = W Exp [(7(0) + V{x)] 



Dynamics 



The arguments above can be readily adapted to the 
exponential potential. By a similar argument as above, 
the dominant contribution to /4l is proportional to if s : 



6 ,(4) 
g J NL 



Oie, 



( x h 



(S 



e e X 
_ '•e'-e ss 
2 £ 2 "c 



0(1) X -V-T?'*. 



(73) 



Once again, in order to damp out non-adiabatic fluctu- 
ations, we need large r/ ss , which in turn implies either 
large 77* or large ?7 X : 



(74) 



The rest of the analysis is identical to the previous section 
and once again leads to exponential damping of of 
the form given in ( 72 ) . 



To change the dynamics one can either look at a wider 
set of trajectories, consider a larger class of potentials, 
or extend the number of inflaton fields. We have already 
discussed the applicability of our results if the inflaton 
motion violates slow roll at the end of section lll Bl While 
we see no reason why similar results to those derived 
above should not hold for more general potentials, we 
have not proved the applicability of our conclusions be- 
yond the class of potentials that we have studied. What's 
more, as detailed in section |IIB[ the particular poten- 
tial we consider is, in some sense, the most general that 
can be readily studied analytically using the SN formal- 
ism. Accordingly, to explore a larger set of potentials we 
would first need to develop alternative (efficient) ways 
of extracting the statistical parameters of interest (the 
power spectrum, the bispectrum etc.). 

For more than two fields, we expect the result to hold 
for appropriately slow roll trajectories in potentials of the 
form wQ>i, ...,<f) n ) = F(Ui(cj) 1 ) + . . . + U n (<j) n )), because 
in this case we can construct a set of conserved quantities 
[24] Ci, . . . ,C„_i by the same method as explained in 
section [TlB| For more general potentials and trajectories, 
no definitive statement can be made. 



The General Case 

A similar argument (though somewhat more involved 
and with a significant loss in clarity as payment for 
the increased generality) holds for the general potential 
W{(f), x) = F(U((j)) + V(x)), and is presented in appendix 

m 



IV. MORE GENERAL MODELS AND OTHER 
APPROACHES TO ADIABATICITY 

Although we have shown that is suppressed dur- 

ing the approach to adiabaticity in a particular class of 
two field inflationary models, it would be helpful to un- 
derstand how generally this conclusion applies. The im- 
portance of this question is made clear when we recall 
that as long as non-adiabatic fluctuations are present, 
the curvature perturbation and its correlation functions 
will evolve outside the horizon. Therefore no sharp pre- 
dictions can be made without an understanding of the 
evolution up to the time of adiabaticity or observation. 
Consequently, no model that relies upon the superhorizon 
evolution of the curvature perturbation to generate non- 
gaussianity can claim to predict a large value for f^i, 
without a calculation which determines its evolution up 
to adiabaticity or observation. 

There are two broad possibilities for generalizing the 
above results: changing the field dynamics or changing 
the mechanism for achieving adiabaticity. 



The Approach to Adiabaticity 

Instead of using a phase of single field inflation to drive 
the cosmological fluctuations into an adiabatic solution, 
one can instead achieve the same result if the universe 
goes through a period of local thermal equilibrium with 
no non-zero conserved quantum numbers (401 I41j . An 
understanding of this process would allow a fuller treat- 
ment of models in which non-adiabatic fluctuations per- 
sist through the end of inflation including more general 
models of multiple field inflation, modulated reheating, 
and the curvaton scenario. We plan to explore this issue 
in future work, and while there is no definitive answer yet, 
we believe that local thermal equilibrium should have a 
similar damping effect as a phase of single field inflation. 

Essentially, this is because the relevant underlying dy- 
namics during a period of local thermal equilibrium are 
very similar to those during a phase of single field infla- 
tion: there is single scalar degree of the freedom (either 
the inflaton or the temperature) which in turn means that 
there are only two solutions for the perturbations in the 
cosmological fluid, both of which must be adiabatic [40] . 
Moreover, the attractiveness of the adiabatic solution as 
one approaches thermal equilibrium is at least qualita- 
tively similar to what happens as one approaches a phase 
of single field inflation. Of course, to show the damping 
(or otherwise) of /nl one needs a model of reheating and 
thermalization to which one can readily apply either an 
analytical method such as the 6N formalism or a suitable 
numerical analysis. 
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CONCLUSIONS 



In this paper, we have considered the generation of 
ft?L from superhorizon evolution of the curvature per- 
turbation in two field inflation, in the most general class 
of potentials (W((j),x) = F{U{4>) + V(x))) that seem 
readily tractable in the 5N formalism. Furthermore, we 
have shown that if the universe passes through a suffi- 
ciently long period of effectively single field inflation to 
force the cosmological perturbations into an adiabatic so- 
lution - as suggested by observation and as is necessary 
to ensure subsequent conservation of the curvature per- 
turbation - I /nl s1 I wn ^ be exponentially damped toward 
a value which is slow roll suppressed. 

We stress, again, that predictions of an observable /nl 
(or any other statistical property of the spectrum) from 
primordial fluctuations are only possible if one knows how 
the fluctuations evolve from the very early universe to 
today. Moreover, due to our ignorance of several phases 
of the early universe, the ability to make sharp predic- 
tions is restricted to the cases where there exists a con- 
servation law which prevents some unknown dynamics 
from changing the statistical properties of the fluctua- 
tions. Luckily, there is such a conservation law when the 
fluctuations are adiabatic. As such, multifield inflation- 
ary models which generate large values for Z^ 3,1 from the 
superhorizon dynamics of the curvature perturbation are 
incomplete without an understanding of the subsequent 
evolution of the spectrum until constancy is achieved. 

To put our results in context, recall that the mea- 
surement of non-gaussianities of the local form has been 
considered to be a sharp probe of the number of infla- 
ton fields as i s small in single field inflation, but 

could apparently grow large in the presence of multiple 
inflatons. This work shows, however, that while it is 
indeed possible to generate a significant contribution to 
the bispectrum in multifield inflation, such a contribution 
is necessarily transitory if one wishes to eliminate non- 
adiabaticity by a phase of single field inflaton. Accord- 
ingly, a significant measurement of /nl 3,1 m the current 
generation of cosmological experiments would present a 
substantial challenge to inflationary model builders. 



Note 



Foundation under Grant No. PHY-0455649 and by the 
Texas Cosmology Center, which is supported by the Col- 
lege of Natural Sciences and the Department of Astron- 
omy at the University of Texas at Austin and the Mc- 
Donald Observatory. J.M. was also supported by the A. 
D. Hutchison Student Endowment Fellowship during the 
completion of this work. 



Appendix A: Details for Homogenous and 
Exponential Potentials 

Homogeneous Potential: W((j),x) = [U(<j>] + V(4>)\ 



With N and C given by (|33j) and (|44j): 

rX 



c = — 2 



N 



p I TT> 

V u 



1 



'Xa 

i r v(x) 



u> (4>) 

Varying N then gives: 
1 



7m 



d X . (44) 



dA = 



1 



™%1 



U'J. \U') C d<f>* \V> 

vj* d x * \u')c °x* \v 



Ax* 



(Al) 



Note that in deriving the above we had to account for 
the dependence of <f> c and \ c on both cj>* and x* ■ We will 
also need: 



dcj) c ( dC 
AC \W* 
dxcfdC 
AC 



dC 

dx* 
dC 

dX* 



Ax, 
Ax* 



From (331 we have: 
dC 1 



dC _ m 2 p 
d~x*~ ~ Vf 



(A2) 



(A3) 



After the first version of this work appeared, a paper by 
Peterson and Tegmark [63] appeared which shows similar 
results, using an alternative framework for analyzing the 
evolution of inflationary perturbations. 
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and 



dC 



" m 



dC _ m 2 p 



dXc V> ■ 

The time t c defines a surface of constant energy: 

W (<j> c , Xc) = constant . 
Differentiating with respect to C then gives: 



A<t> c 
AC 



W 6 



(A4) 



(A5) 



(A6) 



12 



Using the above and W ( p/W x = U'/V, we can differen- 
tiate the expression for C in ( 33 ) and obtain (after some 



manipulation) the following expressions: 



From these it is straightforward to use the relevant SN 
equations to obtain the results given in ( 56 1 and ( 57 ) . 



Exponential Potential: W(<j>,x) = WoExp [U(<f>) + V(x)] 



d(f>c _ 
dC ~ 

AC 



m p 
1 



UL 



V! 



l l 

u? + v? 
1 1 



i u'X 2 

ml U' 2 + V' 2 ' 
m 2 p U' 2 + V' 2 ' 



(A7) 



Substituting (A3) and (A7) into (A2) allows us to read 



off the following: 

c^c = V' 2 U' c V' 2 U' c 

d<p* U' 2 + V' 2 Ui ' d X * U' 2 + V' 2 V.: ' 

dXc = U' c 2 Vj d Xc = U' 2 V' c 
d^* ~ U' 2 + V' 2 Ui ' d X * U' 2 + V' 2 V; 

The derivatives of N are then 



7- (A8) 



dN 1 x h ON 1 y h 



d<f>* m p ' dx* m p ^2eJ 



(A9) 



We have used the slow roll parameters defined in (351; 
while Xh and yh are defined in ( 55 ) : 



1 



Xh 



Vh 



u* + v, 
i 



u* 



Vc4 - U c e* 



Vc4 ~ U c e* 



U*+V* 

In a similar vein, we can find the second derivatives: 
d 2 N 



(55) 



1 



d(j) 2 jm 2 



d 2 N _ 1 

dxl im 2 



d 2 N 1 



rn* 
2d 

,2 



(u c + v c y 1 etc* f 7n , 



n x 



1 - Vh 



2d 



7 + 1 



d(/>*x* 7 m p 



(U c + V c ) 1 efe} ( ir]t 
{U*+V*fd e c 

{Uc + Vc) 2 1 egeg (T& 



- 1 

(A10) 



rj ss is given in ( 54 ) 



(54) 



With N from (46) and C again given by (33), 

rx ! 



C = -~- 2 



1 



''';i. v wm d * + ^J»u%T) dJ(! m 



i 



i 



the variation gives: 
1 



dTV = 



1 



2m 2 p 



2m 2 



1_\ _dh L (±\ _9Xc(l 
U'J 4 \U')c d( t>* \ v 

V'J* d X * \U'Jc d X * \V 



dx* ■ 
(All) 



The analysis begins in the same fashion as in the previous 
section, with the expressions for the derivatives of <f) c and 



Xc with respect to <ft* and x* > given in ( A8 1 . From this 
and ( |A11 1 we obtain the following expressions for the 
first derivatives of N: 



dN 1 x e 



dN 1 y e 



2m p y/M d X* 2m p J^d 



(A12) 



The slow roll parameters are defined in (35). x e and y e 
are defined in (59) as: 



e 4> _ e x 

x e = 1 + 

e c 

c<t> _ f x 



Ue = 1 



The second derivatives are then given by: 



(rf* ~ 24) *e + -^r£ 



(59) 



<9 2 A 


1 1 




d<t>l 


2™ 2 2et 




d 2 N 


1 1 




dx 2 * 




r x 
fc* 


d 2 N 


1 1 




d(j)*x* 


2m2 2ei 





£ 2 'c ' 



4e*e x 

,2 '<= 



(A13) 



?7 SS is given in ( 58 1 : 



»7 



gx^^ _ 4 e e x + e ^x 



(58) 



As above, the SN equations then give the observables of 
J60| and doll. 
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Appendix B: General W(4>,x) = F (U ((f) + V{$)) Results 
1. SN and Observables 

As before we begin with the equations for C and N: 

1 



C = 
N = 



x 1 
, - dy' 

"' xo u'(x') x - 



^IJ X , W x ( X ,<t>(x)) 



dx- 



Thus: 



diV 



2m 2 p 



W 



dcj) c [W_ 
W A 



h 



dXc fW_ 



+ 



2ml 



W 



p i \ " x / * 



W 

dx* 



dx*- (Bl) 



The quantities ii, I2, I3 and I4 are defined as: 

(W x (cf>,x(<t>) 



h 
h 



I A = 



w 4> {<i>,x{4>) 
w^{4>{x)A 



w x (4>{ x ),x 
w x {<t>, x {<t>) 



W*(6x(0) 



W x (cf>(x),X 



w&xWWto (0,xW) \ W) 

Wl^M)) ) dfa 

w 2 (Hx),x)) J i 
W{<f>,x{4>))W^(4>,x{4>))\ 9 x{<t>) 



w 2 (<t>( x ),x)) 



dX: 



dx* 



W(c/ } (x),x)W 4 , x (<f>(x),x)\ d<Kx) 



dx* 



d X . 



(B2) 



Armed with (Bl| and (A8| we can now find derivatives of N. For the first derivatives we have: 



dN 1 


m p 


m p F c Fl4-e* 


8(f)* 2ml 


_V2~4 


V2~d F 'c F * ^ 


dN 1 


m p 




dx* 2ml 







h-h 



h-h 



(B3) 



We've used here (and below) the definitions of the slow roll parameters given in (35 1. For the second derivatives, 
things are a little more complicated: 



d 2 N 1 



dx*d(/)* 



2ml 



d 2 N 1 



dx* 2ml 
d 2 N 1 







Fi 4 


F" 


2d 




\F* 2tt 


'FI 


V* 




(Fi # 


Fi' 


2ei 




{F* 2ei 


Fi 



J> , pip 



€r F'F* 



1 - 



FcFj! 2tt| 
F' 2 + e c 

FcF" , 277- 



pn 



Jx-Jz 



— Je — J& 



2ml 



F*Fi' 



1 ^F c Fif i _F^ L+ 2r ]L] _ j2 _ Ji 



VeUi V F[? J e c \ F* 



(B4) 
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To derive the above we've used the results: 



d<j>* d(f>* 
dx* d\* 



dh 
dl 3 



dh 

d4>* 

dh 



dx* + dx* 



F C F" \ 1 U 2 V' 2 



v* 
YL 



K 2 7 U' 2 U' 2 + V' 2 



+ J1 + J3, 



f*f>j 

F*F" 
F' 2 

F C F" 



1 



1 



F C F"\ 1 U 2 V' 2 



F' 2 
F C F" 
F' c 2 

u' 2 v' 2 



1 



u',v: u' 2 + V 2 

1 u' 2 v' 2 

U'V' U' 2 + V' 2 



+ J2 + J4, 
+ J5 + J7, 



Pi 2 ; V' 2 U' 2 + V' 2 



+ J 6 + h- 



(B5) 



Along with the definitions: 



Ji 



J4 — J7 



+ 



YL 
u> 



+ 



+ 



v 



J 6 = 



-J 



+ 



YL 



h = 



J * 



V' 



'2 



V^f FF^\ V_ 
U> V F' 2 J + U> 



1 - 



FF"\ 9 2 x(0) 



p'2 



d<t>l 



&<t> 



r\ [yll(\- ff " \ 

J* [u> { F' 2 ) ' V 

■2, 



V 



'2 



V f FF"\ d 2 x {4>) 
+ W \ F' 2 ) d<t>*d X * 



d<j) 



'2 



U —(-\- FF "\ — 

v> { "F 2 ) + Y> 



1 - 



FF"\ d 2 (j)(x) 



F '2 



d<j>l 



U' 2 



YLf EEL 
1 ' V 1 f ■■ j 1 r 



,ff^\ d 2 <t>{x) 

F' 2 ) d^d X * 



YL 
u> 



1 



FF' 



+ 



V 2 



F'2 j • JJ, 
FF"\ d 2 X (<py 



p'2 



dx 2 * 



dej) 
U' 2 



ULf FF^ 
V 7 { ^F 2 ) + V 



1 FF"\ d 2 ^(x) 
F' 2 ) d X 2 



dx 



F" FF" 2FF" 2 



F' F' 2 F' 3 



F" FF" 2FF 



"2 



F' F' 2 + F' 3 



F" FF" 2FF" 2 
F 7 F' 2 F' 3 



F FF" 2FF 



"2 



p> F' 2 F' 3 



F FF" 2FF 



p> F' 2 + F' 3 



F" FF" 2FF 

Til ^ 



"2 



F' F' 2 F' 3 



V dfa J 



dx(4>) dx{4>) 
d(j>* dx* 



d(f>* dx* 



(dx(4>)\ 
\ d X * J 



(d<t>{x)\ 
\ d X * ) 



(B6) 



In these expressions functions of x m <A integrals should be rewritten using the solutions x(0) for the particular 
trajectory we are on, so that said integrals are well-defined (and vice- versa in the x integrals). As mentioned in 
the main body of the text, this constrains us to look only at those trajectories where there is a one-to-one mapping 
between (f> and x- 

We can simplify these expressions slightly by defining, in analogy with the homogeneous and exponential cases 
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above 



F F' — f x 
F r F f - e x 



VF = 1 



F'F, e c 



(B7) 



From the above expressions we can (with some work) use the SN formalism to find expressions for P^, and f ] 
in terms of slow roll parameters and derivatives of the potential. First Pq: 



(4) 
NL 



Pr = 



8m4 



(B8) 



Then, for n^: 

tiq — 1 = —2c* 
x M 



1 



"T Zp (7l + 7 2 



1 \ ^<t 



F' 2 



x 



e* e * \ el 



nip 



mp 



(h+h)- 2 + 



\/2eJ 
e* / mp 



(h + h) 



+ 



\ If 9 

x F -x F {h + 7 2 ) 

mp 



Fl 2 2e} 



V2e? 

y F -vf -(h + h) 

mp 



J1+J3+ \I^(J2 + Ji) 



m P 



2eJ 

yp - '— * (h + h) 



mp 



(B9) 



And finally 



6 f (4) 

5 JNL 



ej \ "ip / e * \ m P / 



x 4< 



+ 



Fjy 
F! 2 



f*f" ni 

F? 2ef 



x 2e* 



~x 1 y p 

ei \ mp 



(h + h) 



+ -tt ^f (-fi + h)\ \VF (h + h) 

efe} \ m P \ m P 



5 -^'■i 



KFl_rl_ 
F' 2 2eJ 



-y Uf- ( J 3 + ^4) 



e 
2e? 



\ 2 



(J 6 + Js) |j/f- J ^(/ 3 +/ 4 ) 



(BIO) 



Vdei \ m P I \ m P 



, 4^F c Fj f F C F'J | 2^ 



F' 2 



— hep (7i + 7 2 ) - — \y F (7 3 + 7 4 ) 
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For each of these quantities, if we take the limit t c i* we obtain, as we should, the usual single field results: 



— 1 

6 f (4) 
5 



: /nl 



- -66, + 2r£" , 
= 2e* - rd" . 



(Bll) 



r 



2. The Approach to Adiabaticity and the 
Suppression of /N° L cal 



We will now apply the same analysis as in section III 
to the more general potential discussed above. This case 
is necessarily more complicated than both the homoge- 
neous potential and the exponential potential due to the 
presence of the integrals Ii and Ji as well as terms such 
a s p c ,p* which cannot be expressed in terms of slow roll 
parameters. We will find that these details do not greatly 
affect the conclusions, but the argument is less straight- 
forward in this case. 

We begin by determining the approximate magnitudes 
of the various terms that appear in (BIO I by counting 



factors of slow roll parameters, where as in Section |III| 
e denotes a general first-order slow roll parameter. We 

will take the combinations ' — 



Xf 



-{It +I 2 ) and 



Vf 



(h + h) to be 0(1). This is justified by 



the fact that if either of these combinations were much 
then all of the terms in f-$ would be 



larger, say 0(e 

at most 0(e 2 ), and so would always be small, 
the other hand, if both of these combinations were small 



On 



say 0(e), then the denominator of both (B9| and (BIO) 
would be very small, so in this case, may be large, 
but tiq — 1 will also be large, and so this case is alr ead y 
observationally ruled out. Given the definitions (B7|, 



f,f' 



this also constrains the combination prw- to be < 0(1). 
With these co nside rations in mind, we find that the de- 
nominator of (BIO) is 0(e~ 2 ) 
p39t is 0(e~ r ). 



while the denominator of 



Examining the various terms in ( B9 ) we see that in 



order to maintain a nearly scale-invariant spectrum, we 



must have 



1 



F,F" 



F' 2 



F,F" 



2<i 



F,F" 
F' 2 



27J 



, and 



each of the J j < 0(1). The second through the sixth 
lines of (BIO) are each 0(e _1 ), while the seventh line 
is 0(e~ 2 ). Putting this all together, we see that the 

(4) 

only term in which is not automatically slow roll 



suppressed is the term on the seventh line of (BIO). So 



if /j^l is going to be large, the leading contribution will 
be: 



6 f (4) 

r J NL 



0(e) +0(1) x 



<P f x 



ere 



1 



FcF? 2 V * C 
F' 2 e r 



(B12) 



As discussed in the main text, we wish to examine 
the case where we pass through a phase of effectively 
single field inflation in order to damp away non-adiabatic 
fluctuations. This requires that rj ss becomes large for a 
finite period during inflation, while •q <ja remains small. 
Recalling the definitions ( 35 1 and ( 48 1 we have 



(B13) 



The term - F p, 2 e^e x cannot give the dominant contribu- 
tion to rj ss when ^ ss becomes large, because the same 
term appears in rf a which must remain small otherwise 
inflation will quickly end. Here we have assumed that 
the potential is not so finely-tuned as to create a cancel- 
lation which would keep rf" small while is large. 
As a result, we must have either > \ or r\ x > 
in order to damp away non-adiabatic fluctuations, while 



F '2 



0(e). This means that only the term pro- 
portional to ?7 SS will be important if is large. 

From this point forward, the argument is precisely the 
same as is given in sect ion [ITl] and we find that for rj ss > 1 



Exp 



Hdt 



(B14) 



~ 0(e*) + 0(1) x rTExp 



2 / C v Hri ss dt 



(B15) 



Cjj is a number which is always greater than 1, and de- 
pends on the direction of the effective inflaton during this 
phase (as in section III). Therefore, for any potential of 



the form W(<j>,x) = F{U{4>) + V(x)), we find that // 4) 



NL 



will always be small upon entering the purely adiabatic 
solution after passing through a phase of effectively single 
field inflation. 
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